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We have developed a BFKL-based dipole model in impact parameter space, which has successfully 
predicted a wide range of total cross sections. The model has now been extended to a full event 
generator which is showing promising results. Some of the difficulties in non-leading order BFKL 
are discussed, and preliminary results are compared to CDF and ALICE data. 
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I. INTRODUCTION 

A BFKL-based initial state dipole model in impact pa- 
rameter space has been developed in a series of papers 
[ll-Q. In the Good- Walker approach Q the two incom- 
ing hadronic particles are seen as linear superpositions 
of eigenstates to the interaction. Miettinen and Pumplin 
suggested that the interaction eigenstates are parton cas- 
cades 0|, which in this model are generated with a dipole 
model. Only inclusive, and some semi-inclusive, observ- 
ables have been studied in previous publications, but we 
now want to use the information in the initial state cas- 
cades to produce exclusive final states. Some of the par- 
tons in the initial state cascade will be selected as real 
partons, while others will be seen as virtual fluctuations 
that will be absorbed. The real partons will be further 
evolved through final state radiation with ARIADNE [H 
and hadronisation with PYTHIA H,[l0|]. 

The way the real partons are selected is not straight 
forward, and there are many details in higher orders that 
have large impact on observables. Some of these details 
may be solved by further study, but others seem to be 
out of reach from current perturbative calculations. It 
should be well noted that the results presented here are 
not the only possible way of extending the inclusive for- 
malism. This short summary will not go into detail on 
these complications, leaving them for a publication in the 
near future. 



Here V is the emission probability density, integrated 
over the transverse position z of the emitted gluon. x 
and y are the transverse positions of the partons in the 
emitting dipole. 

The cascade in our model introduces a number of cor- 
rections to Muellers original formulation. 

• The momenta of the partons are necessary to de- 
termine for many of the following corrections. In 
DIPSY, each emission gets a transverse recoil from 
each parent of r t / r\ where r t is the transverse dis- 
tance between the emission and the parent. With 
the transverse momentum set and the rapidity 
known from the emission, the p + of an emission 
can be calculated through 



p+ = p T e v 

where y is the rapidity. Notice that the rapidity is 
negative for particles incoming along the negative 
z-axis. The negative lightcone momentum depends 
on the structure of the rest of the event, but will 
in the cascade be assumed to be that of an on-shell 
gluon: 
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II. THE CASCADE 

The cascades are generated from an initial valence 
state of dipoles in impact parameter space, and is then 
evolved in rapidity with a BFKL-based splitting proba- 
bility. The foundation of the cascade is the leading log- 
arithm formalism first used by Mueller and coworkers 
EMI: 
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To conserve momentum, recoils in px and p + are 
given to the partons in the emitting dipole. The p_ 
will be delivered from the colliding partons from 
the other side, and the limitations are taken into 
account when the interactions are decided. 

Next to leading logarithm effects are known to be 
large in BFKL. While this is not a full NLL calcu- 
lation, it includes what is known to be the domi- 
nating parts of NLL [Ti |. 

— The running coupling can be taken in account 
by replacing a s with a s (/j,) in the emission 
probability. The scale [i is set by the largest 
Pt involved in the emission. 

— The non-singular terms in the splitting func- 
tion suppress large z- values, where the emit- 
ted parton takes the larger part of the momen- 
tum. Most of this effect is included by energy 
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conservation, and requiring that the emitted 
partons are ordered in lightcone momentum 

P+- 

— The energy scale terms are essentially equiv- 
alent to projectile-target symmetry and are 
taken into account by not only requiring or- 
dering in but also in 

• Confinement is added by giving the gluon a mass. 
That modifies the splitting probability to have an 
exponential suppression for large transverse dis- 
tances. This slows down the increase in cross sec- 
tion with energy. 

• Saturation is present in the interaction in Muellers 
original formulation through multiple interactions. 
This only allows for loops that are cut by the inter- 
action frame, not including loops contained in one 
of the cascades. To fully include saturation, and to 
restore frame independence, a 2 to 2 dipole "swing" 
is introduced. The swing can be interpreted either 
as a quadrupole correction, or as the exchange of 
virtual gluons. It tends to swing larger dipoles into 
smaller dipoles, which suppresses the growth of the 
cross section, and can give rise to 2 to 1 mergings 
in the final state, if one of the outgoing dipoles is 
reabsorbed. 

• Coherence is important when a long distance emis- 
sion is made from a parton with other partons close 
by. The long transverse wavelength of the emitted 
gluon cannot resolve the group of closely spaced 
partons, and the emission will treat the group of 
partons as one in terms of recoil and ordering. 
Without coherence, the restrictions of p- ordering 
would be overestimated, and the full phase space of 
allowed emissions would not be taken into account. 



Summing over all pairs of dipoles, one from the state 
incoming from the negative z-axis, and one from the pos- 
itive z-axis, gives the full first order elastic amplitude 
F(b) = J^fij- Allowing multiple interactions gives the 
unitarised elastic amplitude 

T{b) = 1 - e~ F( - b l (3) 

Averaging over the incoming cascades, the interaction 
eigenstates, from both states now gives the elastic cross 
section 

da cl /d 2 b = (* in |T|^ in ) 2 = (^c n 0„|T|^c„(/)„) 2 

= (£ C ^) 2 = < T > 2 (4) 

where is the incoming mass eigenstates and (f> n are 
the interaction eigenstates, that is allpairs of dipole cas- 
cades, with weights c n . T n = 1 — e~ ^*ien ? ij is the scat- 
tering eigenvalue for the pair of cascades n. By replacing 
^in with ^ x the diffractive excitation cross section can 
be calculated, and the optical theorem gives the total 
cross section. From that, the non-diffractive cross sec- 
tion can be found: 

da diB + cl /d 2 b = 5Z(*i„|r|*x>(*x|r|* in ) = (t 2 ). 

X 

da tot /d 2 b = 2(* in |T|* in ) = 2(T) (5) 
do- non _ diS /d 2 b = da tot /d 2 b - da diS+c i/d 2 b = (1 - e~ 2F ), 

where eq. (|3|) was used in the last equality. Identify- 
ing e~ 2F = e~ 2 ^ hi as the non-interaction probability, 
the non-interaction probability of the individual dipole 
pairs factorises, and each pair of dipoles i,j has a non- 
diffractive interaction probability of 1 — e 2 ^ ij . From this 
it is possible to select which dipoles interact in each col- 
lision. 



III. THE INTERACTION 



IV. FINDING THE FINAL STATE 



The elastic scattering amplitude between two collid- 
ing dipoles at a given impact parameter b is in Muellers 
formulation given by 



fij = ./';:i-,.y, ■i- r y j ) = 



log 



(x i -Xj) 2 (y i -y i 



(2) 

This formula is in our model replaced by a similar cal- 
culation in p-space rather than x-space, and modified in 
the same way as the emission in the cascade for recoils 
and lightcone ordering, NLL effects, confinement and co- 
herence. Saturation is already present in Muellers form 
through multiple interactions. The details of these ef- 
fects are mainly guided by frame invariance, requiring 
that the corrections are the same in the interaction as in 
the evolution. 



By tracing the interacting partons back towards the va- 
lence partons, the interacting gluon chains can be found. 
Since the timelike part of the shower will be regarded as 
final state radiation taken care of by ARIADNE, only the 
emissions directly connected to the interacting chains are 
kept. All other emissions are reabsorbed, to avoid double 
counting with the final state radiation. 



A. High pr suppression 



As can be seen from eq. ([IJ dipoles are in the cas- 
cade emitted with a weight of d 2 r/r 2 , corresponding to 
cPpt/Pt- However, the maximum pt in an interaction, 
coming from the cascade or from the interaction, should 
come with a weight d 2 pTjp%- Muellers original model 
was designed for inclusive cross sections, which were not 
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Pseudorapidity distribution at <Js = 1800 GeV 



Charged multiplicity at Vs = 1800 GeV, |/;| < 1, p T > 0.4 GeV 
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FIG. 1: Comparison with a selection of data from CDF [If 
HU. The plots have been generated with Rivet 11711. 




FIG. 2: Comparison with a selection of data from CDF [jgj . 
The plots have been generated with Rivet [l7j |. 



sorbed with a probability 



Pabs = 1 



-PTmin 

2 ' 
^Tmax 



where the maximum px is compared to the closest mini- 
mum along the parton chain. 

The extra phase space opened up by these reabsorb- 
tions (mainly the p^ ordering is significantly changed) is 
approximately covered by the coherence in the cascade: 
A maximum in px means that the reabsorbed dipole is 
smaller than the proceeding emission, and thus that emis- 
sion will not have resolved the small dipole, and the full 
phase space of the merged dipole has been taken into 
account. 



B. FSR and hadronisation 

The gluon chains connected by the exchanges in the in- 
teraction, together with the high-p^ reabsorbtions above, 
determine the ISR. The partons and their colour struc- 
ture are then passed on to ARIADNE [8| to perform 
FSR. The phase space is restricted to spacelike emissions, 
matching the restriction in DIPSY. 

After the FSR, the colour dipoles are allowed to frag- 
ment into hadrons as colour strings using PYTHIA 
[HGIl. Here the low Pt cutoff in ARIADNE is matched 
to PYTHIAs starting scale to cover all phase space with- 
out double counting. 



RESULTS 



affected by small dipoles, and this discrepancy was not 
a problem. To produce a correct p^-spectrum in the fi- 
nal state though, it is necessary to reweight the maxima. 
Reweighting according to 

[18], 

maxima in px are reab- 



As can be seen in previous publications, the following 
inclusive cross sections are well described at high energy: 

• Total and elastic pp as function of \fs and t, 
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FIG. 3: Comparison with a selection of data from ALICE 




• Single and double diffractive excitation in pp as 
function of Mx , 

• Total 7*p as function of W and Q 2 (also low Q 2 
with VMD), 

• Deeply virtual compton scattering as function of 
W, Q 2 and t (also low Q 2 with VMD), 

• j*p — > pp as function of W, Q 2 and t (also low Q 2 
with VMD). 

The exclusive final state extension introduces many am- 
biguities from higher orders or non-perturbative QCD, 
which cannot easily be solved from first principles. All 
the plots are from a preliminary version of DIPSY from 
spring 2010. We will here discuss how it compares to 
data, and how sensitive the observables are to details in 
the algorithm. 

A selection of data is shown in figs [TJ [5] and [31 

• The pseudorapidity distribution fits well with data 
in this preliminary version of DIPSY, but it is sen- 
sitive to details regarding mainly how the p± or- 
dering is implemented. 

• The leading px distribution in region close to the 
trigger particle fits well with data and is relatively 
stable between different versions. A successful pre- 
diction. 

• The distribution in A(f> around a minijet does not 
have enough activity in the transverse and away 
region. This is related to the rapidity distance be- 
tween a minijet and its recoil. Possibly it is caused 
by an overestimate of the absorbtion of partons due 



to local pT-rnaxima, giving too large rapidity dis- 
tances between the remaining partons. 

• The multiplicity distribution from CDF follows 
data approximately over several orders of magni- 
tude. This holds true for most version of DIPSY. 
The bias towards low multiplicities in this version 
is possibly related to the same overestimate of re- 
absorbtion. The same effect is seen in the ALICE 
multiplicity distributions, but at the higher ener- 
gies DIPSY undershoots data somewhat in the tail. 

The suppression algorithm for high px for this prelimi- 
nary version of DIPSY has a known flaw, allowing some 
small dipoles to survive unsupressed. This gives a too 
strong tail in the px spectrum, but is corrected in later 
versions. 

dNch/dr) in midrapidity at ALICE overshoots simula- 
tions at high energies in this version, but it is sensitive 
to several details in the model. More analysis is needed 
to understand the energy dependence in DIPSY. 



VI. CONCLUSIONS AND OUTLOOK 

We have over the last years developed a BFKL based 
dipole model in transverse impact parameter space, in- 
cluding the major parts of NLL, confinement and satu- 
ration. This model has in previous publications shown 
to reproduce a wide range of total, elastic and diffractive 
cross sections with few parameters. Now we are produc- 
ing exclusive final states with the same model. There 
are however many effects in higher orders, which it is not 
clear how to include, but that still has a significant effect 
on observables. Some effects can be accounted for with 
intuitive solutions, while others require tuning from data. 
The model can reproduce many exclusive observables at 
the Tevatron and the LHC. It should be noted that the 
results here are from a preliminary version of DIPSY, 
with several known flaws, and we expect the results in 
an upcoming publication to be improved. 

While all data here are from pp colliders, the model can 
be applied to any high energy collisions, as for example 
DIS, AA, pA and eA. AA will be the first to be studied 
after pp. 

We would in the future also like to extend the model 
to include diffractive final states, but since the interac- 
tion amplitude does not factorise as neatly as in the non- 
diffractive case, it is more complicated. We do have some 
ideas though, and hope to return to this topic in future 
publications. 
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